Galactic Point Sources of TeV Antineutrinos 
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High energy cosmic ray experiments have identified an excess from the region of the Galactic 
Plane in a limited energy range around 10^* eV (EeV). This is very suggestive of neutrons as 
candidate primaries, because the directional signal requires relatively-stable neutral primaries, and 
time-dilated neutrons can reach Earth from typical Galactic distances when the neutron energy 
exceeds an EeV. We here point out that if the Galactic messengers are neutrons, then those with 
energies below an EeV will decay in flight, providing a flux of cosmic antineutrinos above a TeV 
which is observable at a kilometer-scale neutrino observatory. The expected event rate per year 
above 1 TeV in a detector such as IceCube, for example, is 20 antineutrino showers (all flavors) and 
a 1° directional signal of 4 events. A measurement of this flux can serve to identify the first 
extraterrestrial point source of TeV antineutrinos. 

PACS numbers: 13.85.Tp, 95.85.Ry, 95.30.Cq 



An intriguing anisotropy in the cosmic ray spectrum 
has emerged in the energy range near an EeV. The Akeno 
Giant Air Shower Array (AGASA) has revealed a corre- 
lation of the arrival direction of the cosmic rays to the 
Galactic Plane (GP) at the 4(T level l]. The GP excess, 
which is roughly 4% of the diffuse flux, is mostly con- 
centrated in the direction of the Cygnus region, with a 
second spot towards the Galactic Center (GC) Evi- 
dence at the 3.2a level for GP enhancement in a similar 
energy range has also been reported by the Fly's Eye Col- 
laboration . The existence of a point-like excess in the 
direction of the GC has been confirmed via independent 
analysis 'i^ of data collected with the Sydney University 
Giant Airshower Recorder (SUGAR). This is a remark- 
able level of agreement among experiment using a variety 
of techniques. 

Independent evidence may be emerging for a cosmic ac- 
celerator in the Cygnus spiral arm. The HEGRA exper- 
iment has detected an extended TeV 7-ray source in the 
Cygnus region with no clear counterpart and a spectrum 
not easily accommodated with synchrotron radiation by 
electrons @. The difficulty to accommodate the spec- 
trum by conventional electromagnetic mechanisms has 
been exacerbated by the failure of CHANDRA and VLA 
to detect X-rays or radiowaves signaling acceleration of 
any electrons p. The model proposed is that of a pro- 
ton beam, accelerated by a nearby mini-quasar or pos- 
sibly Cygnus X-3, interacting with a molecular cloud to 
produce pions that are the source of the gamma rays. 
Especially intriguing is the possible association of this 
source with Cygnus-0B2, a cluster of more than 2700 
(identified) young, hot stars with a total mass of ~ 10'' 
solar masses ■ Proton acceleration to explain the TeV 
photon signal requires only 0.1% efficiency for the con- 
version of the energy in the stellar wind into cosmic ray 
acceleration. Also, the stars in Cygnus-0B2 could be the 



origin of time-correlated, clustered supernova remnants 
forming a source of cosmic rays. By cooperative acceler- 
ation their energies may even exceed the ~lPeV cutoff 
of individual remnants and accommodate cosmic rays up 
to the ankle, where the steeply faUing (oc i<;-3-i6±o.08^ 
cosmic ray spectrum flattens to £'-2.8±o.3 

All evidence points to a transition from galactic to ex- 
tragalactic sources above several EeV of primary energy. 
The steepness of the falloff between the knee (about 3 
PeV) and the ankle (about 10 EeV) is expected from 
supernova shock models, and may indicate that we are 
witnessing the high energy end of the galactic flux. The 
extension of the nominal PcV cutoff beyond the ankle can 
be understood as a collective effect of stellar winds orig- 
inating in the region of multiple supernova explosions. 
These provide a second acceleration to the particles and 
boost their energies far beyond the values expected from 
their single-shock encounter. Strong additional support 
for this picture emerges from accumulating evidence (in 
a bi-modal proton-iron model) for a dominant Fe compo- 
nent in the flux 9] . The importance of the heavy compo- 
nent is apparent all the way down to the region of several 
PeV 1^, with the spectral index hardening slightly to 
3.02 ± 0.03 below 500 PeV. An immediate consequence 
of this nucleus- dominance picture is the creation of free 
neutrons via nuclei photodisintegration on background 
photon fields. These liberated neutrons are presumably 
responsible for the observed directional signals. This im- 
plies that it may not be a coincidence that the signal 
appears first at energies where the neutron lifetime al- 
lows propagation distances of galactic scales, i.e., 10 kpc. 

For every surviving neutron at ~ EeV, there are 
many neutrons at lower energy that decay via n — > 
p + e~ -\-Ve. The decay mfp of a neutron is cr„r„ = 
10 (Sn/EeV) kpc, the lifetime being boosted from its 
rest-frame value r„ — 886 seconds to its lab value via 
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r„ = En/rUn- The proton is bent by the Galactic mag- 
netic field, the electron quickly loses energy via syn- 
chrotron radiation, and the V^, travels along the initial 
neutron direction, producing a directed TeV energy beam 
whose flux is calculable. We show in this Letter that 
the expected v Hux from the direction of the Cygnus 
region is measurable in IceCube ^Jl- Furthermore, by 
the same logic, a km-scale Mediterranean detector, if 
designed with sufficiently low threshold, can see the D 
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The variables appearing in Eq. |^ are the antineutrino 
and neutron energies in the lab (EV and En), the antineu- 
trino angle with respect to the direction of the neutron 
momentum, in the neutron rest-frame (9—), and the an- 
tineutrino energy in the neutron rest-frame (e-) . The last 
three variables are not observed by a laboratory neutrino- 
detector, and so are integrated over. The observable E-jy 
is held fixed. The delta-function relates the neutrino en- 
ergy in the lab to the three integration variables. The 
parameters appearing in Eq. arc the neutron mass 
and rest-frame lifetime (m„ and t„), and the distance to 
the neutron source (_D). dFn/dEn is the neutron flux at 
the source, or equivalently, the neutron flux that would 
be observed from the source region in the absence of neu- 
tron decay. Finally, (e— ) is the normalized probability 
that the decaying neutron in its rest-frame produces a 
T7g with energy e-p-; de-p- = 1 defines the nor- 

malization, where the maximum neutrino energy in the 
neutron rest frame is just Q = mn—mp—rrie = 0.71 MeV, 
and the minimum neutrino energy is zero in the massless 
limit. ""^ For the decay of unpolarized neutrons, there is 
no angular dependence in 

The En integration in Eq. is effectively cut off at 
~ EeV, the energy beyond which a neutron is stable over 
a 10 kpc path-length, and it is truly cut off by the i?max 
of the neutron spectrum. The expression in parentheses 
in Eq. is the decay probability for a neutron with 
lab energy En, traveling a distance D. In principle, one 
should consider a source distribution, and integrate over 
the volume J d'^D. Instead, we will take D to be the 1.7 
kpc distance from Earth to Cygnus 0B2; for the purpose 
of generating the associated neutrino flux, this cannot be 



The massless-neutrino approximation seems justifiable here: 
even an eV-mass neutrino produced at rest in the neutron rest- 
frame would have a lab energy of nii, Fn ^ GeV, below threshold 
for neutrino telescopes. 



flux pointing toward the GC source as well. The GZK 
neutrinos |l2l | and the "essentially guaranteed" D flux cal- 
culated here probably constitute the best motivated cos- 
mic neutrino fluxes. Of these two neutrino fluxes, the 
expected event rate for the galactic beam is higher: 4 9^ 
events per year and 16 in 4- f'r showers. 

We turn to the calculation. The basic formula that 
relates the neutron flux at the source {dFn/dEn) to the 
antineutrino flux observed at Earth {dFir/dEv) is: 



dP , , d cos 077 ^ ^ / — N , n 

d^ (^-) j ^ ^ - ^" (1 + M/^n] . (1) 
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in error by too much. 

The Galactic anisotropy observed by the various col- 
laborations spans the energy range 0.8 to 2.0 EeV.^ The 
lower cutoff specifies that only neutrons with EeV en- 
ergies and above have a boosted ct„ sufficiently large to 
serve as Galactic messengers. The upper cutoff reflects an 
important feature of photodisintegration at the source: 
heavy nuclei with energies in the vicinity of the ankle 
will fragment to neutrons with energies about an order 
of magnitude smaller. To account for the largest neutron 
energies, it may be necessary to populate the heavier nu- 
cleus spectrum in the region above the ankle. This is 
not a problem - one fully expects the emerging harder 
extragalactic spectrum to overtake and hide the steeply 
falling galactic population. It is not therefore surprising 
that in order to fit the spectrum in the anisotropy region 
and maintain continuity to the ankle region without in- 
troducing a cutoff, the AGASA Collaboration required a 
spectrum cx E~^ or steeper 1]. 

A detailed scenario for ultrahigh energy nuclei (par- 
ents of the anisotropy neutrons) originating in a pulsar 
close to the Cygnus 0B2 region has been recently de- 
scribed Q. The sequence begins with the one-shot ac- 
celeration in the spinning neutron star |l5j . resulting in 
an E~^ spectrum. Softening of the spectrum to 
ensues through gravitational wave losses during spin- 
down [l^ . Following this scheme, we assume that some 
of the nuclei are captured in the dense region of the 
source, attaining sufhcient diffusion in milliCauss mag- 
netic fields. The resulting time delay of several thou- 



^ Actually, the anisotropy reported in 3] peaks in the energy bin 
0.4 — 1.0 EeV, but persists with statistical significance to energies 
as low as 0.2 EeV. The full Fly's Eye data include a dire ctio nal 
signal which was somewhat lost in unsuccessful attempts Il3l to 
relate it to 7-ray emission from Cygnus X-3. This implies that 
if neutrons are the carriers of the anisotropy, there needs to be 
some contribution from at least one source closer than 3-4 kpc. 
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sand years[l7| produces a further steepening of the injec- 
tion power law spectrum. Note that once diffusion has 
been established, additional Rayleigh steps in the Galac- 
tic magnetic field do not change the spectral index signif- 
icantly. In their random traversal of the OB association, 
the nuclei undergo photodisintegration on far infrared 
thermal photons populating molecular clouds with tem- 
peratures of 15-100 K y^]. Taking an average photodis- 
integration cross section of 40 mb, we find an interaction 
time between 4 and 1300 yr, allowing sufficient neutron 
production to explain the anisotropy. 

To incorporate the preceding discussion in our work, 
we take in what follows a single power law neutron 
spectrum with a spectral index of 3.1, representing an 
average over the PeV-EeV energy region. Specifically, 
dFn/dEn = CE~^-^, with the normalization constant 
fixed near an EeV to the observed excess. The constant C 
is determined by integrating dFn/dEn over a bin {Ei , E2) 
with the result 

Fn= C E-^-^ = 0.95 C En^'^ sinh [1.1 A] , (2) 

J El 



where En = VWE^ ~ lO^'^ GeV, A = \n{E2/Ei) w 
1.38, and F„ « 9 km'^ yr"! j^. 

The typical energy for the antineutrino in the lab is 
that of the decaying neutron times Q/m„ ^ lO^'^. Thus, 
the stability of neutrons at ^ EeV implies a PeV up- 
per limit for the produced antineutrinos. The increasing 
abundance of neutrons below an EeV in turn implies an 
increasing neutrino flux as energies move below a PeV. 

Nuclei with Lorentz factor ^ 10^ are synthesized in all 
supernovae. Hadronic interactions with the HII popula- 
tion (density 30 cm""^ and photodisintegration from 
ultraviolet photons emitted from OB stars results in a 
flux of PeV neutrons. From the measured Lyman emis- 
sion of the O stars 0, the 0-^B temperature, and lu- 
minosity characteristics of the B stars we obtain a 
photon number density in the ultraviolet of ^ 230 cm^"^ 
for the core region ~ 10 pc. This implies a collision 
time of about 0.1 Myr, comparable to the time scale for 
hadronic interactions. However, in contrast to hadronic 
interactions, significant photodisintegration occurs in the 
outer 20 pc of the source. The photon density is reduced 
to > 25 cm~'^, which lengthens the reaction time to ~ 
1 Myr. The diffusion time 1.2 M yr) is a bit smaller 
than the age of the cluster ^2.5 Myr [20j, and somewhat 
higher than the reaction time. This is sufficient to per- 
mit over 90% efficiency for photodisintegration over the 
lifetime of the source. The effective volume for photodis- 
integration is then about a factor of 27 larger than for 
hadronic interactions, and the net result of all these con- 
siderations is that the PeV neutron population is about 
an order of magnitude greater than TeV charged pions 
resulting from hadronic collisions . 

Performing the cos 6-jj -integration in Eq. |Q over the 



delta-function constraint leads to 

dEn dFn 



dF- nin f 
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En dEn 
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(£;„)(l_e-a) 
(3) 



^min _ Errrrin „„J pn 



An approxi- 
mate answer is available. Setting the beta-decay neu- 
trino energy e— equal to its mean value = ep, we have 
^(cf) = S{e—~ eo).^ When the delta-function is substi- 
tuted into Eq. one gets 



dEn dFn 



^ En dE, 



{En){l-e 



dE;^"^' = 2^0 

(4) 

Further approximation is available. Treating the neu- 
tron decay factor, 1 — exp(---), as a step function 
Q^^max _ ^^^^ energy E^'"'' - 0{^^) = 

(10^) EeV, i.e. the neutron is unstable for En < E™^'^ 
and stable for En > E™'^^, one obtains from Eq. Q 
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We have found that for E'^'''' = 3 EeV, Eqs. ^ and Q 
differ by less than 1% for Eu as high as 1 PeV. 

A direct Ve event in IceCube will make a showering 
event, which, even if seen, provides little angular resolu- 
tion. In the energy region below 1 PeV, IceCube will re- 
solve directionality only for and P^. Fortunately, neu- 
trino oscillations rescue the signal. Since the distance 
to the Cygnus region greatly exceeds the Vf. oscillation 



length Ao 



10 



-2 Eb 
PoV 



parsecs (taking the solar oscilla- 



tion scale Sm? ^ 10~^eV^), the antineutrinos decohere in 
transit. The arriving antineutrinos are distributed over 
flavors, with the muon antineutrino flux i^p^ given by the 
factor i sin^(2 00) ~ 0.20 times the original F^^ flux. 
The Vr flux is the same, and the flux is 0.6 times the 
original flux. Here we have utilized for the solar mix- 
ing angle the most recent SNO result 9q ~ 32.5° ji^ . 
along with maximal mixing for atmospheric v^^-Vr neu- 
trinos and a negligible z/e component in the third neutrino 
eigenstate. 

The integral neutrino flux Fy{> E-) = ^^_dEjT^i 
is particularly useful for experiments having a neutrino 



^ The delta-function in the neutron frame gives rise to a fiat spec- 
trum for the neutrino energy in the lab for fixed neutron lab- 
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energy _E„ = r„ m„ : 

dP _ dcosH- f dt-\ ( dP" 

dE^~J-i 2 \dE^) \d^J 2r„eo 

with < ET7<2r„to. 
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FIG. 1: Integrated flux of (solid) and j/p + J^e + t'r (dashed- 
dotted) predicted to arrive at Earth from the direction of the 
Cygnus region. Also shown is the integrated Vf^ + V/j. atmo- 
spheric flux for an angular bin of 1° x 1°. The shaded band in- 
dicates the region excluded by the AMANDA experiment . 
The expected number of showers i/p (solid) and + Vs + Vt 
(dashed-dotted) to be detected (say in IceCube) are plotted 
on the bottom-left. The expected background for the same 
angular bin is indicated by the dashed line. 



detection-efficiency that is independent of neutrino en- 
ergy, or nearly so. IceCube is an example of such an 
experiment. Our calculated integral flux is shown in 
Fig.n As mentioned above, the nuclear photodisintegra- 
tion threshold leads to an infrared cutoff on the primary 
neutron energy at the source, which in turn leads to a 
low energy cutoff ~ TeV on the integral flux. 

We now estimate the signal-to-noise ratio at IceCube. 
The angular resolution of the experiment « 0.7° allows a 
search window of 1° x 1° 0|. We begin with the "noise" . 
The event rate of the atmospheric z/-background that will 
be detected in the search bin (Afiioxio ~ 3 x 10^'' sr) is 
given by 



dN 



background 



Aeff y" dEJ^+i,{E) p{E) An 



l°xl° , 



(6) 

where A^s is the effective area of the detector, J,j^ij{E) 
is the -\- atmospheric flux in the direction of the 
Cygnus region (about 40° below the horizon) and 
p{E) « 1.3 X 10^6 (£'/TeV)°-8 denotes the probabihty 
(generic to ice/ water detectors) that a v (or v) with en- 
ergy E on & trajectory through the detector produces 
a signal [i^. For a year of running at IceCube and 
^min _ Y TeV, from Eq. lO one obtains a background 
of 1.5 events. Existing limits on the 7-ray flux from the 



Cygnus region j23| provide an upper limit on neutrino 
fluxes generated via tt^ decay at the source. This limit 
is below the atmospheric background in the region of in- 
terest. Poisson statistics then imply that a signal > 3.7 
events is significant at the 95% CL |2^. The number of 
V showers in the signal, for energies above i?™™, is given 
by 



dN 
~dt 



signal 



/■ dF- 

Acff / dE^—f{E^)p{E^)^ 

7£;mi„ dET7 



(7) 



For a year of running at IceCube, one expects 20 neutrino 
showers (all flavors) with energies > 1 TeV, of which 4 
events will cluster within 1° of the source direction, 
comfortably above the stated CL. Neutrino flux at 1 TeV 
may also originate in the decay of 1 PeV neutrons from 
sources whose spectrum cuts off at that energy, and hence 
are not subject to normalization by the anisotropy. Thus 
our estimate may be regarded as very conservative. 

IceCube is not sensitive to TeV neutrinos from the 
GC, as these are above the IceCube horizon, where atmo- 
spheric muons will dominate over any signal. However, 
other kilometer-scale neutrino detectors, such as those 
planned for the Mediterranean Sea, may see the GC flux.* 
Additionally, Southern Auger should see the cosmic ray 
excess in the direction of the GC (29|] and Northern Auger 
should be sensitive to the Cygnus region. 

We conclude that in a few years of observation, Ice- 
Cube will attain 5cr sensitivity for discovery of the Fe— *■ 
71 — !■ i7g — !■ cosmic beam, providing the "smoking ice" 
for the CP neutron hypothesis. 
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